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ABSTRACT 

Recently, Park and Gott reported an interesting observation: image separation of 
lensed QSOs declines with QSO redshift more precipitously than expected in any 
realistic world model, if the lenses are taken to be either singular isothermal spheres 
or point masses. In this Letter I propose that the observed trend arises naturally if the 
lensing galaxies have logarithmic surface mass density profiles that gradually change 
with radius. If the observed lack of central (odd) images is also taken into account, 
the data favor a universal dark matter density profile over an isothermal sphere with 
a core. Since the trend of image separation vs. source redshift is mostly a reflection of 
galaxy properties, it cannot be straightforwardly used as a test of cosmological models. 
Furthermore, the current upper limits on the cosmological constant may have to be 
revised. 

Key words: gravitational lensing - cosmology - galaxies: structure 



1 INTRODUCTION 

Since gravitational lensing takes place over cosmological dis- 
tances, it can in principle be used to measure cosmologi- 
cal parameters, such as mass density, Q, and cosmological 
constant, A. For example, the frequency of multiply imaged 
QSOs is a sensitive function of A, and therefore the observed 
abundance of lensed QSOs has been used to place upper lim- 
its on A (Turner 1990; Fukugita et al. 1992; Kochanek 1996). 
Another example is the relation between image separation, 
8im and source redshift, z 3 of multiply split QSOs, which 
mostly depends on the sum of f2 and A, and thus should 
serve as a good indicator of the curvature of the Universe. 
(Gott, Park & Lee 1989; Fukugita et al. 1992, Park & Gott 
1997, hereafter PG97). If the Universe is flat and galaxies 
are singular isothermal spheres, the di m -z„ relation should 
be flat, i.e. independent of source redshift. If the Universe 
is open, there should be a mild decline of image separations 
with redshift; a ~ 20% decline is expected between z s ~ 1.5 
and 5, in an extreme fi = case. 

The mildness of the expected 9i m -z s trend can be ex- 
plained as follows. For any given lens, and fixed lens and 
source redshifts the image separation of images with compa- 
rable brightness is quite independent of source impact pa- 
rameter. Since optimal lens redshift for any high redshift 
source is roughly the same, zi ~ 0.5, and angular diame- 
ter distances vary little past z ~ 0.5, the observed image 
separation should change little with z s . If a population of 



lenses is considered, having a range of properties, like scale 
lengths and velocity dispersions, there will be a correspond- 
ing spread in image separations, but still little or no trend 
with z s - 

Yet the observed 9i m of lensed QSOs declines strongly 
with z s (See Figure 1 of PG97, or Figure || of this paper.) 
PG97 present a list of source redshifts and image separations 
of 20 multiply imaged QSOs. The source redshifts span a 
range between 1.4 and 4.5, and image separation between 
0.5" and 7". There is a drop of at least a factor of 4 in image 
separation from the row-z s to high-z s cases. According to 
PG97, in a flat Universe such distribution can be ruled out 
at 99% confidence level, while in the most extreme empty 
Universe, it is ruled out at 97%. To make the confidence 
level decrease below 95%, one would need to 'remove' 2-3 
largest separation or highest redshift cases from the sample. 
(See PG97 for a discussion of statistical significance, possible 
errors, etc.) 

In this Letter, I suggest that the trend, which is due 
to the lack of high redshift wide separation lensed QSOs is 
naturally reproduced if a set of three conditions is met, (i) 
galaxies' central mass profiles have logarithmic slopes that 
change with radius, (ii) there is a dispersion in galaxy prop- 
erties, like central surface mass densities or velocity disper- 
sion, and (iii) the characteristic length scale of galaxy dark 
matter halos varies with galaxy luminosity as r oc L a , where 
< a<0.5. 
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2 MODEL 

In this section I describe the lensing model, and introduce 
some simplifying assumptions about the lenses, galaxies, and 
sources, QSOs. I do not aim to produce a detailed model be- 
cause the small number of lensed QSOs in the current sample 
does not necessitate it. I need not assume any particular cos- 
mological model, because the effect implied by the present 
model is considerably stronger than that of cosmology. 



2.1 Model Assumptions 

2.1.1 Lenses 

I consider two types of galaxy mass profile, an isothermal 
sphere with a core, ISC, and a universal dark matter pro- 
file, NFW, derived from numerical simulations of Navarro, 
Frenk & White (1995, 1996). These profiles were chosen be- 
cause their logarithmic surface mass density gradually flat- 
tens with decreasing radius. The reason for using these, 
as opposed to constant power law profiles, like a singular 
isothermal sphere, is discussed in Section H. Both ISC and 
NFW provide a plausible description of the dark matter ha- 
los of galaxies. ISC is commonly used because of the ob- 
served shapes of the rotation curves of spiral galaxies. Addi- 
tionally, based on the results of the HST Snapshot Survey, 
Maoz & Rix (1993, hereafter MR93) conclude that ellipticals 
must possess ISC-type dark matter halos in order to repro- 
duce the observed image separation distribution of multiply 
imaged QSOs. The NFW profile has been recently claimed 
as a universal dark matter halo profile describing halos of 
gravitationally bound structures over 4 orders of magnitude 
in mass. It is also supported on theoretical grounds (Evans & 
Collett 1997). I assume that all galaxy- lenses are circularly 
symmetric. 

A description of the lensing properties of these profiles 
can be found in Schneider et al. (1992), Bartelmann (1996), 
and Williams & Lewis (1997). It suffices to say here that 
ISC is roughly fiat within core radius r c , and has an isother- 
mal density profile outside. The vertical normalization of 
the profile is fixed by its central surface mass density, kq, 
in terms of critical density for lensing, E cr ; t . NFW is sin- 
gular at the centre, its logarithmic slope gradually steepens 
from p oc r _1 at the centre, to p oc r -3 at large radii, and 
is isothermal around the characteristic scale length, r s . The 
vertical normalization of its projected surface mass density 
is proportional to k s , which is given in terms of E CJ -jt- 

As was pointed out by PG97, the scatter in the observed 
Oim-Zs points is large. One of the main sources of scatter is 
probably the spread in galaxy-lens properties. For example, 
a singular isothermal sphere lens of a given a v produces a 
constant bending angle at the lens, and hence image sepa- 
ration varies as Di s /D oa , where the D's are the lens-source 
and observer-source angular diameter distances. If all lens- 
ing galaxies had the same a v , one would expect the observed 
angular image separations to be proportional to the corre- 
sponding Di 3 /D oa values, whereas a plot of Di 3 /D os vs. 8i m , 
in 8 systems where both z s and Zi are known reveals no cor- 
relation at all, and a scatter in 9i m of a factor of a few larger 
than the scatter in Di s /D oa . This strongly suggests that the 
spread in galaxy-lens properties is important. Therefore I as- 
sume a family of galaxy-lenses, each having the same scale 
length, r c and r s for ISC and NFW respectively, but a range 



of ko and k s , which is equivalent to a range in velocity dis- 
persions. 

I assume that r c and r s and galaxy luminosity function 
do not evolve with redshift, at least in the interval most 
relevant here, i.e. optimal lens redshift for high-z sources, 
zi ~ 0.3 — 1. This assumption is supported by the recent 
spectroscopic observations of Lilly et al. (1995), which indi- 
cate that the luminosity function of elliptical galaxies, i.e. 
the galaxy population that is believed to provide the bulk 
of the lensing optical depth (see MR93), evolves very little 
from 2 ~1 to the present. 

When the source impact parameter is sufficiently small, 
all the NFW lenses, and ISC lenses with ko > 1 produce 
3 images, whereas the observed multiply imaged systems 
mostly have an even number of images, either 2 or 4. This 
is thought to be because the central (odd) image is demag- 
nified below visibility. I will return to the importance of 
the central odd image later, in Section 2.2. The circularly 
symmetric lenses considered here cannot produce 5 image 
systems; observed 4(+l) cases are a property of elliptical 
lenses. However the present treatment will not suffer if only 
symmetric lenses are considered. 

For sources at large redshifts, the optimal lens redshift 
increases very slowly with z s in all cosmologies except for 
the most extreme A dominated cases. So I will assume that 
zi is constant, independent of z s . The critical surface mass 
density for lensing, E cr it changes little with z s , if zi is fixed, 
and sources are at high z s . For example, in an £7=0.3 open 
Universe, moving the source from z s =l to 5 reduces E cr i t by 
17%, if zi = 0.2. Therefore, I will assume that galaxy-lens 
projected surface mass density is given by ko and k s , with 
no dependence on redshifts. 

2.1.2 Sources 

Analytic fits to QSO luminosity function (QLF) usually take 
the form of a double power law, with steep bright-end slope 
and shallow faint-end slope; the transition occurs at some 
characteristic luminosity, Lo, which evolves with redshift, 
and is commonly parameterized as Lo(z s ) oc (l + z 3 ) a . Boyle 
et al. (1990) derive a pure luminosity evolution for z<2 
QSOs, with a=3.2, while Hewett et al. (1993) conclude that 
QLF changes shape with redshift, and that evolution must 
slow down beyond z ~ 1.5 compared to the predictions of 
Boyle et al. For 2<z<3 they derive a ~1.5. QSO evolu- 
tion at higher redshifts, up to 5 is less constrained, though 
it is sometimes assumed that the shape of the QLF does 
not evolve beyond z ~ 2, while Lo evolves such that it 
'slides' brightward along the high luminosity part of the 
QLF (Wallington & Narayan 1993). I assume that at any 
given redshift there are no QSOs brighter than Lo, and that 
QLF at L < L is N QS o(L) oc L' s , with s = 1.2 (MR93). 
I adopt the a parameterization of QSO evolution for high 
redshifts, and derive results for a=0, and 2. 



2.2 Results 

We are interested in two lensing properties of the galaxy- 
lenses: image separation and total image magnification. 
Given a galaxy-lens, a source can have a range of impact 
parameters, each resulting in a different image magnifica- 
tion, and image separation, even though the latter tends to 
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stay rather constant, roughly equal to twice the Einstein 
ring radius of the lens. To account for a range of source im- 
pact parameters, I calculate the impact parameter weighted 
averages of image magnification, (fi) and image separation, 
(8) as a function of ko or k s , for ISC and NFW respectively. 

Figure |l| shows the results. Vertical normalization of im- 
age separation is irrelevant for now. The important feature 
is that for both profiles the image separation increases with 
ko or k s , while total image magnification decreases. In other 
words, given a population of galaxy-lenses with a range of Ko 
or k s , image magnification goes inversely with image sepa- 
ration. Now, image magnification couples to source redshift, 
as follows. Flux of an Lo QSO is given by, 

/ oc n(z a )Lo(z s )/D L (z a ) 2 

oc n(z s )(l + z s ) a /D A {z 3 ) 2 {l + z a ) 4 

oc [i(z a )(l + z s ) a ~ 4 , (1) 

where j-i(z s ) is total image magnification, and Dl and Da 
are the luminosity and angular diameter distances respec- 
tively. In the second and third lines I used the fact that 
Dl — (1 + z) 2 Da oc (1 + z) 2 , since angular diameter dis- 
tance stays roughly constant for high-z s sources. So to make 
it into a magnitude limited sample, with a certain fu m , an 
Lo QSO at z s has to be magnified by at least 

fl(zs) OC flim(l + Z a ) 4 ~ a , (2) 

with fainter QSOs magnified more. Because larger redshift 
QSOs undergo larger magnifications, a fixed fu m of any 
given survey translates into a lower limit on (/i) as a func- 
tion of z s , which in turn implies an upper limit on {#), as 
read off from Figure |l|. The resulting relation is plotted as 
solid (ISC) and dashed (NFW) curves in Figure g, for two 
values of QSO luminosity function evolution parameters, a; 
0, and 2. Both sets of curves reproduce the upper envelope 
of the observed points quite well. The paucity of observed 
points does not allow to differentiate between a of and 2, 
or between ISC and NFW profiles. 

The curves are scaled vertically and horizontally to 
roughly match the observed distribution. Adjusting the scal- 
ing in the vertical direction yields the angular size of r c 
and r 3 respectively. For ISC angular core radius is 0.37", 
while for NFW angular scale radius is 1.25". For typical lens 
redshifts of 0.5 these translate into r c = 1.4/i _1 kpc, and 
r s = 4.6/i _1 kpc. These values are roughly what one would 
derive from modelling individual multiple- image systems us- 
ing ISC/NFW profiles. Scaling the curves in the horizontal 
direction gives the average minimum magnification of lensed 
QSOs as a function of z 3 . For a — models, {fx) / fu m = 
0.039(1 + z s ) 4 , and for a = 2, (fi) / fum = 0.4(1 + z s ) 2 ; i.e. if 
QSO characteristic luminosity does not evolve, a lower limit 
on magnification of QSOs at z = 3 is (fJ,)/fu m =W, whereas 
it is 6.4 if Lq(z s ) oc (1 + z 3 ) 2 . Surveys with brighter fu m have 
higher (fi), but the dependence on (l + z s ) is weak, and most 
existing surveys have fu m of 1-2 magnitudes within each 
other. 

Let us summarize the conclusions thus far. Roughly, 
Figure ^ implies that average image magnification is in- 
versely proportional to the galaxy central surface mass den- 
sity, (fA) oc k,q].; and that average image separation is di- 
rectly proportional to the galaxy central surface mass den- 
sity, {0) oc Ko,s- In general, the observed image separation 
is also proportional to the characteristic scale length of the 
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Figure 1. Magnification and image separation vs. the central 
surface density of two different types of lenses: ISC and NFW. 
Both magnification and image separation arc impact parame- 
ter weighted averages of these quantities. The dotted line is the 
impact parameter weighted average magnification of a singular 
isothermal sphere. Note that for ISC and NFW profiles the mag- 
nification and image separation go in the opposite sense with 
increasing Ko and k s , while for SIS the average magnification is 
always a constant (see Section ^) . 

galaxy-lens, so (8) oc Ko, 3 r c ,s, but since I have so far assumed 

that all galaxies have the same characteristic scale length, I 

have not used the latter dependence. The scale length can 

be a function of central concentration, r c oc Kq s . Combin- 

_ i 

ing these relations gives, (fj) oc (8) 1 +f , and together with 
eq.(2) it implies, 

(0)oc(l + z s )- (1+/3)(4 ~ a) (3). 

With (3 = this equation basically reproduces the curves 
in Figure || Note that (3 does not have to be 0, any value 
between roughly -0.5 and 2 will do, i.e. the total exponent in 
eq. (3) should be around -few. I will return to the discussion 
of j3 and its implications for the dark matter halos of galaxies 
in Section ^| 

Thus both ISC and NFW galaxy-lens models can re- 
produce the observed #i m -z s trend, and neither is preferred 
based on these observations alone. However, another widely 
documented observation can break the tie; all multiply im- 
aged QSOs lack the central odd-numbered image, i.e. ob- 
served cases are either doubles or quadruples. This has long 
been interpreted as evidence in favor of centrally condensed 
galaxy profiles, which would demagnify the central image 
below visibility. The other possible explanation, dust obscu- 
ration in the centres of galaxies, is largely ruled out because 
radio observations of multiply split QSOs also do not reveal 
central images (Myers et al. 1995). I apply this argument 
here: ISC model predicts that kq ~ 1 — 2 models should 
have central images of brightnesses comparable to the pri- 
mary image, while with the NFW models the central images 
of low K a lenses should be 5-10 times demagnified compared 
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Figure 2. Image separation vs. source redshift. The points are the 
observed multiply imaged QSOs (taken from Table 1 of PG97); 
stars represent the lenses from the HST Snapshot Survey (MR93), 
circles are from other, ground based surveys. Predictions of the 
isothermal sphere with a core, ISC and universal dark matter 
halo, NFW models (labeled) wer e no rmalized vertically and hori- 
zontally, as described in Section 2.2. Both ISC and NFW profiles 



can account for the observed trend. 



to the primary image. This rules strongly in favor of the 
NFW model. 



3 IMPLICATIONS FOR GALAXY MASS 
PROFILES 

The key feature of the galaxy lens population that allows 
to reproduce the envelope is that image separation goes in- 
versely with image magnification. This allows magnification 
bias, which acts through image magnification, to couple im- 
age separation to source redshift. This property of a lens 
population is a direct consequence of the individual galax- 
ies having a changing logarithmic slope, as in the ISC and 
NFW models. Lenses that do not have this property, for ex- 
ample singular isothermal spheres and point masses, cannot 
reproduce the observed envelope. 

For a singular isothermal sphere (SIS) the average im- 
pact parameter weighted magnification can be calculated to 
be 4, independent of the velocity dispersion, a v , while image 
separation is proportional to a 2 . So (fi), being a 'universal' 
constant, does not couple to (8), and thus SIS cannot repro- 
duce the observed trend. The property that (fj,) is constant 
is not unique to SIS, but can be generalized to any family of 
mass profiles with no in-built scale length, like single power- 
laws and point masses. In such models the only parameter is 
the absolute mass normalization of the profile, which deter- 
mines the size of Einstein ring radius, and hence the image 
separation. However, (fj.) in such self-similar profiles is a con- 
stant, independent of absolute mass normalization. 

Therefore galaxies must have changing logarithmic sur- 



face mass density profiles in order to be able to reproduce 
the di m -Zs relation. 



4 DISCUSSION 

There are two interesting aspects in the observed distribu- 
tion of image separations vs. source redshift. In this Letter I 
set out to explain one of these, the lack of high redshift wide 
separation lenses, i.e. the Qi m -z s anticorrelation. It turns out 
that the model presented in Section [| also naturally explains 
the other interesting feature of the 9i m -z a plot, namely the 
existence of wide separation cases; there are 8 lensed QSOs 
with 9im > 3". Most models currently found in the literature 
have trouble predicting a large population of wide separation 
lenses; for example MR93 predict the peak in 6i m distribu- 
tion at < 1", with vanishingly small number of cases above 
3". So it has been argued that large 8i m cases are not due to 
isolated galaxies, but are the result of cluster-aided galaxy 
lensing. QSO 0957+561, with 6.1" between its two images, 
is adduced as supporting evidence. However, if cluster-aided 
galaxy lensing is the correct explanation, then wide sepa- 
ration lenses should be found at all redshifts, and not just 
at low z s , as is currently the case (Figure ^). Furthermore, 
PG97 show that with the help of cluster lensing, image sep- 
aration should increase with z s . Since this is clearly not ob- 
served, cluster aided lensing is probably not important in 
general (0957 must be a special case) , and the large number 
of wide separation lenses needs an explanation. 

Let us derive the distribution of image separations as 
predicted by the model of Section ^. (I will use the ISC 
model because all the derivations required in this section 
can be carried out analytically with this model. However 
the general arguments presented here will also apply to the 
NFW profile.) So far only the dark matter halo properties 
have been discussed. These need to be related to the observ- 
able galaxy properties, like luminosity. I will adopt relations 
similar to those used in MR93 (their Section 2.2.1); 



oc L a , M/L oc L b , hence kq oc L 1 ' 



(4). 



Optical effective radii of ellipticals are observed to scale as 
L 1 ' 2 (Lauer 1985), and since it is commonly assumed that 
they are linearly proportional to the galaxies' dark mat- 
ter scale length, a is usually taken to be 1.2. From the 
observations of the fundamental plane of ellipticals (Kor- 
mendy & Djorgovski 1989), b w 0.25. Similar assumptions 
were made by Kochanek (1996). Parameter (3 of eq. (3) is 

= a/(l - 2a + &|J 

The frequency of multiply imaged QSOs as a function 
of image separation is given by, 



dF 

W) 



dN(L) 
dL 



cVr(L)}- 



dL 



d(6(L)) { fum J 



(5) 



where dN(L)/dL oc L~ 12 e~ L is the Schechter luminosity 
function (Schechter 1976); L is in units of L*, and I assume 
a constant slope of -1.2. The lensing cross section for a galaxy 
of luminosity L is given by the radial caustic in the source 



t Note that MR93 and Kochanek (1996) assumed a=1.2, and 
b = 0.25, and hence effectively their j3 = —1.04, which would not 
reproduce the 9i m -z s anticorrelation, see eq. (3). 
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plane, whose radius is y r (L) = (k^ 3 — 1) 3//2 . The image 
separation is approximately equal to the diameter of the 
tangential critical curve in the lens plane, (6) ~ 2r c \/ Kq — 1. 
Using eq. (4), 

d(e(L)} 2r c ,,L a -\Kl[l-a + b]-a) 

dL («g - 1) ' 1 j 

Here, r Cj * and ko,* refer to the core radius and kq of an L* 
galaxy. I take r c> * to be 1.4/i _1 kpc, equal to the constant 
core radius derived in Section The value of Ko,* is then 
calculated from r c ,* and an assumed asymptotic line of sight 
velocity dispersion, a, — 300km s _1 ; no,* = 5. The last term 
in eq. (5) is the magnification bias. 

The function dF/d(8) of eq. (5) is plotted in Figure |[ 
for a range of a corresponding to the allowable range of /3, in 
eq. (3), o=0.0, 0.2, 0.4, and 0.6; b = 0.25 was used through- 
out. The angular splitting by an L* galaxy is denoted by an 
arrow. These predicted distributions apply to the source red- 
shift range where the upper envelope, described by eq. (3) 
and plotted in Figure ^ does not cut in, i.e. for z B < 2.5. The 
solid histogram is the observed distribution for the same 
source redshift cutoff. Even though the small number of cur- 
rently observed lenses does not allow to make any precise 
conclusions, it is apparent that a ~ 0.4 reproduces the his- 
togram quite well. The corresponding value of /3 is 0.9, which 
is perfectly consistent with the allowed range of /3, -0.5 to 
2 (see eq. [3]). Thus, a ~ 0.4 (j3 ~ 0.9,) reproduces both the 
dim-Zs anticorr elation, and the observed frequency of wide 
separation lenses. In this scenario, the 9i m ~ 7" lenses are 
due to ~ 2L* galaxies, with a mass within 10/i _1 kpc of 
3.7 1O 11 /!," 1 M . 

The major deviation of the present model from those 
found in the literature is in the value of a, the power law 
index relating the dark matter characteristic scale length of 
a galaxy to its luminosity. As mentioned earlier, a is usually 
assumed to be 1.2, while consistency with lensing observa- 
tions in the framework of the present model implies a ~ 0.4. 



5 CONCLUSIONS 

If the observed distribution of image separation vs. source 
redshift of lensed QSOs is real and not a result of false 
lenses, like physical QSO pairs, then it has the following im- 
plications for galaxy dark matter halos: (i) galaxies must 
have changing logarithmic surface mass density profiles, 
with a 'universal dark matter' model being preferred over an 
isothermal sphere with a core, (ii) there must be a spread in 
galaxy properties, like the central surface mass density, (iii) 
the characteristic length scale of dark matter halos should 
scale with luminosity as L a , where a ~ 0.4. 

The corollary of the present model is that higher z B mul- 
tiply imaged QSOs are statistically more magnified, and are 
preferentially lensed by galaxies with lower central surface 
mass densities and intrinsic luminosities. These predictions 
can be tested if each lens case is modelled individually, which 
will become possible when the lensing galaxies of most mul- 
tiply split QSOs are detected. 

The strong dependence of the observed 6 im -z s relation 
on galaxy properties means that it cannot be easily used as 
a test for the curvature of the Universe unless we know a lot 
more about the mass distribution of galaxies. Furthermore, 
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Figure 3. Predicted image separation distribution for lensed 
QSOs. The curves are labeled by parameter a, defined by r c oc L a . 
Vertical normalization is arbitrary. The predictions are valid for 
z s where the magnification bias effect, discussed in Section 
and shown in Figure ^, does not impose an upper limit on (?; m of 
high-z s lensed QSOs, i.e. for z s <, 2.5. Accordingly, the observed 
distribution to which the curves apply, is restricted to z s < 2.5 
and is shown as the solid histogram. Values of a ~ 0.4 seem to 
describe the observed distribution adequately. In these models 
image separation is a monotonically increasing function of galaxy 
luminosity; image separation produced by an L* galaxy is indi- 
cated by an arrow. 

because of the implications the present model has on the 
lensing properties of galaxy population, like the dependence 
of the predicted frequency of lensed QSOs on parameter a 
(see Figure ^), the current upper limits on A may have to 
be revised. 
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